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In order to study if the bulk viscosity may induce a big rip singularity on the flat FRW cosmologies, 
we investigate dissipative processes in the universe within the framework of the standard Eckart the- 
ory of relativistic irreversible thermodynamics, and in the full causal IsraelStewart-Hiscock theory. 
We have found cosmological solutions which exhibit, under certain constraints, a big rip singularity. 
We show that the negative pressure generated by the bulk viscosity cannot avoid that the dark 
energy of the universe to be phantom energy. 
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I. INTRODUCTION 

The existence of a exotic cosmic fluid with negative 
pressure, which constitutes about the 70 per cent of the 
total energy of the universe, has been perhaps the most 
surprising discovery made in cosmology. This dark en- 
ergy is supported by the astrophysical data obtained from 
Wilkinson Microwave Anisotropy Probe (WMAP) (Map) 
and high rcdshift surveys of supernovae. 

The dark energy is considered a fluid characterized by 
a negative pressure and usually represented by the equa- 
tion of state w = p/ p, where w lies very close to —1, most 
probably being below — 1 . Dark energy with w < — 1 , the 
phantom component of the universe, leads to uncommon 
cosmological scenarios as it was pointed out in . First 
of all, there is a violation of the dominant energy condi- 
tion (DEC), since p + p < Q. The energy density grows 
up to infinity in a finite time, which leads to a big rip, 
characterized by a scale factor blowing up in this finite 
time. These sudden future singularities are, neverthe- 
less, not necessarily produced by a fluids violating DEC. 
Barrow has shown, with explicit examples, that exist 
solutions which develop a big rip singularity at a finite 
time even if the matter fields satisfy the strong-energy 
conditions p > Q and p -\- ip > 0. A generalization of 
Barrow's model has been realized in Q, giving its La- 
grangian description in terms of scalar tensor theory. It 
was also proved by Chimento et al ^ that exist a duality 
between phantom and flat Friedmann-Robertson- Walker 
(FRW) cosmologies with nonexotic fluids. This duality is 
a form-invariance transformation which can be used for 
constructing phantom cosmologies from standard scalar 
field universes. Cosmological solutions for phantom mat- 
ter which violates the weak energy condition were found 
in i. 



The role of the dissipative processes in the evolution of 
the early universe also has been extensively studied. In 
the case of isotropic and homogeneous cosmologies, any 
dissipation process in a FRW cosmology is scalar, and 
therefore may be modelled as a bulk viscosity within a 
thermodynamical approach. 

A well known result of the FRW cosmological solu- 
tions, corresponding to universes filled with perfect fluid 
and bulk viscous stresses, is the possibility of violating 
DEC 0. The bulk viscosity introduces dissipation by 
only redefining the effective pressure, P.J^^ , according to 

p,, =p + n = p-3eH, (1) 

where 11 is the bulk viscous pressure, ^ is the bulk vis- 
cosity coefficient and H is the Hubble parameter. Since 
the equation of energy balance is 



p + 'iH{p + p + TT) = 0. 



(2) 
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the violation of DEC, i.e., p + p + Ii < implies an in- 
creasing energy density of the fluid that fills the universe, 
for a positive bulk viscosity coefficient. The condition 
^ > guaranties a positive entropy production and, in 
consequence, no violation of the second law of the ther- 
modynamics |3|- 

In the present paper we show that the above results are 
straightforward to obtain from the exact cosmological so- 
lutions already found by Barrow in 6]. These solutions 
were obtained using non causal thermodynamics. Never- 
theless, we consider a more physical approach like the full 
Israel-Stewart-Hiscock causal thermodynamics, showing 
that it is also possible to obtain big rip type solutions. 

The organization of the paper is as follows: in Section 
II we present the field equations for a flat FRW universe 
filled with a bulk viscous fluid within the framework of 
the Eckart theory. We indicate that under a constraint 
for the parameters of the fluid, one of the Barrow's so- 
lutions presents a future singularity in a finite time. In 
Section III we obtain big rip solutions using the approach 
of the full Israel-Stewart-Hiscock causal thermodynam- 
ics. In Section IV we discuss our results in relation to 
the nature of the dark energy of the universe. 
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II. ECKART THEORY 

The FRW metric for an homogeneous and isotropic flat 
universe is given by 



'dt^ + a{if {dr^ + r^{de^ + sin^0d(j>^)) , (3) 



where a{t) is the scale factor and t represents the cosmic 
time. In the following we use the units SttG = 1. In the 
first order thermodynamic theory of Eckart ig] the field 
equations in the presence of bulk viscous stresses are 



P 
3' 



with 



and 



a 



n = -3H^. 

The conservation equation is 

p + 3H{p + p + n) = 0. 



(4) 
(5) 

(6) 
(7) 
(8) 



Assuming that the dark component obey the state equa- 
tion 



P= il- 



(9) 



where < 7 < 2, we can obtain from equations (0)) to 
a single evolution equation for H: 



2H + = 3^H. 



(10) 



This equation may be integrated directly as a function 
of the bulk viscosity. For 7 9^ the solution has the form 



Hit) = 



62 



imdt 



C + §7 /ei 



(11) 



where C is an integration constant. From this equation 
we find the following expression for the scale factor 

2/(37) 



from which we conclude that p = 3iJ^ + const. Compar- 
ing this expression with Q we have that the integration 
constant is zero. 

Thus we have that for the state equation p — — p, i.e. 
for 7 = 0, the scale factor is not defined by the field 
equations. So for a given a{t) we can write H and then 
obtain the expressions for the energy density from Eq. Q 
and the bulk viscosity from Eq. H13I) . Clearly if ^ = the 
well known de Sitter scale factor a{t) — e^o* is obtained, 
where p = —p and both are constants. 

Notice that the solution of the field equations may be 
written through ^(i) or a{t) because there are three inde- 
pendent equations for the four unknown functions a{t), 
p{t), C(t) andp(i). 

Now we are interested in the possibility that there are 
cosmological models with viscous matter which present 
in its development a big rip singularity. 



A. The case for 7 7^ 

Firstly, let us consider the case 7^0. If the viscous 
fluid satisfies DEC, then the condition < 7 < 2 must be 
satisfied. Thus for 7 < we have a phantom cosmology. 
Now from the thermodynamics we know that ^ > 0, and 
if 7 < the Eq. 1)12(1 implies that we can have a big rip 
singularity at a finite value of cosmic time. 

Let us consider some examples to see this more clearly. 
From Eq. H12|) the well known standard case for a per- 
fect fluid, i.e. ^ = 0, takes the form a{t) = D{C + 
(3/2)7f)2/(3T). This scale factor may be rewritten as 



a{t) 



1 



2/(37) 



and the energy density is given by 

Pa 



P = 



;i+i^o7t)^ 



(15) 



(16) 



where p^ — 3H^^, in order to have H{tg — 0) — H„ > 0. 
If 7 < we have a big rip singularity at a finite value of 
cosmic time t^^^ — —2/{3Hgj) > =0. 

In the special case of ^{t) — — const we have 
from Eq. (fT^ for the scale factor a{t) = D{C + 
(t/Co) e'^^/^)^o*)2/(37)_ \Ye can rewrite it into the form 



ait) = D[c+-^J e'^f ^W"^* dtj , (12) a{t) = a„ (^1 + |l 7 (e ^ */2 _ 1) j 



2/(37) 



where _D is a new integration constant. Thus for a given 
^{t) we have the expressions for a{t), p{t) and p{t). 
For the case 7 = we have from Eq. HlOfl 

and substituting this expression into Eq. ||SJ) we have 

p = 6HH, (14) 



from which we obtain for the energy density 



pit) = Po 



(l + :|i^(e3?o */2_l)y 



(17) 



(18) 



where = 3 H^. As before, for 7 < we have a big rip 
singularity at a finite value of cosmic time t^^^ = ln{l~ 

7fc)>*o=0. 
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Note that any additional condition on the system of 
the field equations will fix the unknown functions. So 
for instance, for a variable ^{t) we can take the condition 

Another example in this line is given by the solution 
obtained by Barrow ^ for the case ^ ^ p^l"^ . Effectively, 
Barrow @ assumed that the viscosity has a power-law 
dependence upon the density 



a > 0. 



(19) 



where a and s are constant parameters, and exact cos- 
mological solutions for a variety of ^(p) in the form given 
by equation ((T^ . In particular, for the case s — 1/2, 
i.e., ^ — ap^^"^ , yields a power-law expansion for the scale 
factor. Nevertheless, none condition was imposed upon 
the parameters a and 7 in order to obtain solutions with 
big rip. 

For the case s — 1/2, the integration of equation IjlOfl 
yields the following expression for H{t) 



1 



1 



(Via - 7) - h) : 



(20) 



where Hq — H(t — to) and to correspond to the time 
where dark component begins to become dominant. The 
scale factor becomes 



a{t) = flo 1 



tbr 



(21) 



where ao — a{t = to)- If we demand to have the occur- 
rence of a big rip in the future cosmic time then we have 
the following constraint on the parameters a and 7 



VSa > 7, 



(22) 



leading the scale factor blow up to infinity at a finite time 
tbr > to, which expression is 



(23) 



In terms of time ti,r, the Hubble parameter is given by 



Hit) = i/o 1 



t-tp 



(24) 



From the equation I^J and the parameterized equations 
(|21|1 and H24|) for the scale factor and Hubble parameter, 
respectively, we obtain the expression for the increasing 
density of the dark component in terms of scale factor 



p(a) = 



3(V3a-7) 



(25) 



We reproduce completely this solution if we put into the 
Eq. p2fl the bulk viscosity given by 



^{t)^VSaHo(l- 



t-to 

^br 



B. The case for 7 = 



(26) 



Notice that the structure of equation l|l()|l changes if 
7 = and ^ is an arbitrary function of the density, since 
the quadratic term in H dissapears. Nevertheless, if ^ '--^ 
p^/^, the structure of the equation ((Tn|l is the same for 
any value of 7 in the range < 7 < 2, except in the case 
VSa = 7, where equation 1)1 0|l becomes H = 0. Then, 
the solution with 7 = can be obtained directly from 
the general solution given by equation (|21|l . In this case 
there is a big rip singularity at a finite value of cosmic 
time i,^ = 2/{3V3H„a) > t^ = 0. 



III. ISRAEL-STEWART-HISCOCK THEORY 

We now consider the dissipative process in the uni- 
verse within the framework of the full causal theory of 
Israel-Stewart-Hiscock. In this case we have the same 
Friedmann equations but instead of equation Q, we have 
an equation for the causal evolution of the bulk viscous 
pressure, which is given by 



rtl + Il = -3^H-Kul3H+^-^-^y (27) 



where T is the temperature and r the relaxation time. 
In order to close the system we have to give the equation 
specifying T 

T = pp''. (28) 
The the relaxation time is defined by the expression 

T=^=ap'-\ (29) 
P 

where /3 > 0. This model imposes the constraint 



1 



7 



(30) 



in order to have the entropy as an state function. Notice 
that the above constraint exclude the range < 7 < 1, 
which implies that quintessence fluids are not allowed 
in this approach. With the above assumptions the field 
equations and the causal evolution equation for the bulk 
viscosity lead to the following evolution equation for H 
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if+-(l + (l-r)7) HH 



(1 + r) H-^H^ + ^ (7 - 2)) i/^ + i32-^a-i7i/ 



4-2s 



0. (31) 



r 



As in the non causal case we will choose s — \/2 and the 
above equation becomes 

H + bHH - (^2--^ H-'^H^ + aH^ ^0, (32) 



where a is defined by 



V3c 



and b by 



6=3 1 



V3c 



(33) 



(34) 



Solutions of equation l|32|l .were obtained in In this 
work only was considered 7 in the range 1 < 7 < 2. Some 
of these solutions presents an increasing energy density 
and accelerated expansion. 

Inspired in the solution for the Hubble parameter given 
by equation H24() in the non causal scheme, we use the 
following Ansatz, where for simplicity we take tg = 



H{t) = A{nr~t) \ 



(35) 



where A = HoTbr. With this Ansatz the scale factor a{t) 
evolutes as 



a{t) {ur - ty 



(36) 



and the energy density, p, of the dark component as a 
function of the scale factor becomes 



pia) ~ a^/^. 



(37) 



Using the Ansatz H35|) in equation H32|l we obtain a 
second grade equation for A 

aA'^+bA+- = 0. (38) 
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The solutions for A are given by 

2A± = --±VA, (39) 
a 

where the discriminant A has the expression: 

Since we are interested only in positive solutions for A, 
the coefficient a must be negative. We have two cases of 
interest. 



Case a < 0; 7 > 0. In this case only A+ correspond to 
a solution with big rip. The parameters a and 7 satisfy 
the following constraint 



VSa > 7 (1 - I) 



(41) 



Notice that there is no big rip solution if the cosmic fiuid 
representing the dark component is stiff matter (7 = 2). 
The factor (l — f ) is the correction introduced by the 
causal thermodynamics to the constraint given by equa- 
tion (|22ll . The solution for is given by 



A+ = 



1 



V3q 



^ + 2 



1/2 



1- 1 + 



' V3ay 2 

which implies that a big rip will occurs at a time 



Ur = A+Hr. 



(42) 



(43) 



The expressions for a — a(t) and p = p{a) can be easily 
evaluate from equations H36|) and (|37|l . respectively. 

Case a < 0; 7 < . Since we need A > in order to 
have real solutions, the parameters a and 7 must satisfy 
the following constraint 



^/3o< 



(44) 



If A = 0, i.e., VSa = ^\ 7 | /2, the solution for A, 
which we shall call Aq, has the following expression 



\/3q 



If A > 0, the solutions for A can be written as 
A± = Ao ± i\/A. 



(45) 



(46) 



IV. DISCUSSION 

Within the framework of the non causal thermody- 
namics we have showed that the power law solution, 
found by Barrow in 6] for dissipative universes with 
^ = ap^/'^, yields cosmologies which present big rip sin- 
gularity when the constraint given in equation H22|l holds. 
If we consider that the dark component is quintessence, 
i.e., < 7 < 2/3, with a sufficiently large bulk viscosity 
will make this quintessence behaves like a phantom en- 
ergy. In the range 2/3 > 7 > 2 it is possible, at least 
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from the mathematical point of view, to obtain sohitions 
with big rip even with a matter fluid. It is not clear for 
us how can be interpreted a radiation fluid, for exam- 
ple, with a large bulk viscosity leading to high negative 
pressures and increasing densities. 

At the boundary between the quintessence sector and 
the phantom sector, i.e. 7 = or p = — p, also there exist 
cosmologies with a big rip singularity. 

Using a more accurate approach like the the full causal 
theory of Israel-Stewart-Hiscock, we have also found cos- 
mological solutions with big rip. 

If 1 < 7 < 2 the parameters a and 7 satisfy the con- 
straint given in equation H41() . Due to the constraint 
stiff matter is not allowed. As we mentioned above this 
correspond to matter fluids that can lead to a phantom 
behavior. Quintessence region are not allowed. If 7 < 
the cosmological solutions can be computed directly from 
equations (|45(l and (|46|l . The main conclusion, in the 
context of the full causal thermodynamics, is that in or- 
der to obtain physically reasonable big rip solutions, the 
dark component must be phantom energy. 

Note added: While this manuscript was being written 
we noticed about the work of Brevik and Gorbunova . 
The authors also consider the possibility of big rip in vis- 
cous fluids with p — wp, by a different formalism. They 
consider the case where the bulk viscosity is proportional 
to the scalar expansion. This is equivalent to the Bar- 



row's choice ^(t) cx -^p, and then they also reobtained 
the Barrow's solution, for 7 7^ 0, considered here. In the 
framework of the standard Eckart theory citeEckart, the 
authors show that fluids which lie in the quintessence re- 
gion (w > —1) can reduce its thermodynamical pressure 
and cross the barrier w = — 1 , and behave like a phantom 
fluid {w < —1) with the inclusion of a sufficiently large 
bulk viscosity. The case for 7 = was not considered by 
these authors. 
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